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Abstract 

An  experimental  study  to  determine  the  through-plane  thermal  conductivity  of  dry  Nafion®,  various  diffusion  media,  catalyst  layer,  and  the 
thermal  contact  resistance  between  diffusion  media  and  a  metal  plate  as  a  function  of  temperature  and  pressure  was  performed.  Dry  Nafion® 
thermal  conductivity  was  determined  to  be  0.16  ±  0.03  W  m_1  K_1  at  room  temperature,  which  decreases  to  0.13  ±  0.02  Wm~*  K_1  at  65  °C. 
Diffusion  media  thermal  conductivity  was  found  to  be  function  of  PTFE  content  and  manufacturer,  and  was  0.48  ±  0.09  W  m  1  K_1  for  untreated 
and  0.22  ±  0.04  W  m-1  K_1  for  20  wt.%  PTFE  treated  SIGRACET®  diffusion  media,  respectively.  Toray  diffusion  media  thermal  conductivity  was 
measured  to  be  1.80  ±  0.27  Wm~‘  K_1  at26°C  and  decreases  to  1.24  ±  0.19  Wm-1  Kat  73  °C.  The  thermal  contact  resistance  between  Toray  carbon 
paper  and  aluminium  bronze  material  was  determined  to  vary  from  6.7  x  10-4  to  2.0  x  10~4  m2  KW_I  for  an  increase  in  compression  pressure  from 
0.4  to  2.2  MPa.  The  equivalent  thermal  conductivity  of  a  0.5  mg  cm~2  platinum  loaded  catalyst  layer  was  estimated  to  be  0.27  ±  0.05  W  m~*  K_1 .  A 
one-dimensional  analytical  model  was  also  used  to  estimate  the  temperature  drop  in  the  fuel  cell  components.  A  maximum  of  3-4  °C  temperature 
drop  can  be  expected  for  a  200  p,m  thick  SIGRACET®  diffusion  media  at  1  A  cm~2.  The  thermal  properties  characterized  should  be  useful  to  help 
modelers  accurately  predict  the  temperature  distribution  in  a  fuel  cell. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Detailed  knowledge  of  the  internal  temperature  distribution 
in  a  polymer  electrolyte  fuel  cell  (PEFC)  is  critical  for  efficient 
water  and  thermal  management.  In  a  fuel  cell,  local  variation 
in  temperature  can  be  attributed  to  the  waste  heat  generation, 
which  includes  the  irreversible  heat  of  electrochemical  reaction, 
losses  due  to  over-potential  at  each  electrode  and  Joule  heating 
in  all  components,  but  mostly  the  electrolyte.  The  temperature 
distribution  can  affect  the  drying/flooding  and  degradation  phe¬ 
nomenon  in  the  fuel  cell,  which  deteriorates  its  performance. 
Therefore,  a  thorough  knowledge  of  in  situ  temperature  dis¬ 
tribution  is  essential  for  understanding  the  thermal  and  water 
transport  in  a  PEFC. 

Direct  in  situ  temperature  measurement  introduces  additional 
challenges  due  to  the  minute  length  scale  involved,  anisotropic 
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nature  of  porous  media,  and  existence  of  complex  two-phase 
flow  in  the  PEFC.  Various  modeling  approaches  [1-4]  have 
been  presented  to  predict  the  temperature  distribution  of  PEFC 
under  different  operating  conditions.  However,  limited  infor¬ 
mation  available  in  the  literature  on  the  varying  through-plane 
and  interfacial  thermal  resistance  values  of  fuel  cell  components 
presents  a  barrier  for  the  modelers.  Highly  accurate  estimation  of 
temperature  distribution  is  impossible  without  complete  knowl¬ 
edge  of  the  cell  component’s  thermal  properties  under  different 
operating  conditions,  PTFE  content,  and  compression. 

Several  studies  have  been  performed  by  various  researchers 
to  determine  the  thermal  conductivity  of  various  materials  for 
non  fuel  cell  applications.  Washo  and  Hansen  [5]  used  the 
steady  state  method  to  measure  thermal  conductivity  of  oriented, 
amorphous  specimens  of  polymethylmethacrylate  (PMMA)  and 
polystyrene.  In  the  field  of  micro-electro-mechanical  system 
(MEMS),  numerous  studies  have  been  performed  to  deter¬ 
mine  the  thermal  property  of  thin  films  and  polymers  [6-15]. 
Kurabayashi  et  al.  [8]  have  used  steady  state  Joule  heating 
and  harmonic  Joule  heating  method  with  electrical  thermom- 
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etry  resistance  (or  optical  thermometry  [9]),  to  measure  the 
in-plane  and  through-plane  thermal  conductivity  of  spin-coated 
polyimide  films.  Choy  et  al.  [10,11,16]  have  developed  a  flash 
radiometry  technique  to  determine  the  thermal  conductivity 
(thermal  diffusivity)  of  polymer  films.  Kurabayashi  [6]  have  pre¬ 
sented  a  brief  review  of  various  experimental  and  mathematical 
models  for  measurement  of  thermal  conductivity  in  the  in-plane 
and  through-plane  directions  for  various  polymers.  A  guarded 
heat  flow  meter  (using  ASTM  E1530  procedure)  was  used  by 
Song  et  al.  [17,18]  to  measure  thermal  conductivity  of  lithium 
polymer  electrolyte  and  composite  cathode  for  lithium  batter¬ 
ies.  Morotta  and  co-workers  [19,20]  also  studied  the  thermal 
conductivity  and  thermal  contact  resistance  of  various  thermo¬ 
plastic  and  thermosetting  polymers  using  the  flux  measurement. 
A  detailed  review  is  presented  by  Savija  et  al.  [21]  on  various 
analytical  and  experimental  methods  for  thermal  contact  resis¬ 
tance. 

To  the  author’s  knowledge,  none  of  the  above  methods  for 
thin  films  have  been  used  for  determination  of  thermal  conduc¬ 
tivity  of  fuel  cell  materials.  Vie  and  Kjelstrup  [22]  estimated 
thermal  conductivity  of  a  fuel  cell  membrane  by  measuring  in 
situ  temperature  gradients  in  a  fuel  cell.  Burford  and  Mench 
[23]  also  estimated  thermal  conductivity  of  diffusion  media 
and  Nafion®  by  measuring  electrolyte  temperature  using  micro¬ 
thermocouples  embedded  in  the  electrolyte.  Non-isothermal  and 
two-phase  modeling  literature  [1-3,24,25]  has  used  inconsistent 
thermal  conductivity  values  for  PEFC  materials.  Table  1  sum¬ 
marizes  the  thermal  conductivity  values  of  fuel  cell  material 


Table  1 

Thermal  conductivity  (W  m  1  K1 )  of  fuel  cell  components  reported  by  various 
authors  in  the  literature 


Diffusion  media 

Catalyst 

layer 

Membrane 

Author/reference 

1 

0.5,  1.0,2.94“ 

- 

0.95a 

Ju  et  al.  [1] 

2 

- 

- 

0.67a 

Beming  et  al.  [2] 

3 

1.7“ 

- 

0.5a 

Hwang  [3] 

4 

- 

- 

0.1±0.1b 

Vie  and  Kjelstrup  [22] 

5 

0.13— 0.19b 

- 

0.025-0.25b 

Burford  and  Mench  [23] 

6 

1.5a 

1.6a 

0.34a 

Rowe  and  Li  [24] 

7 

0.15a 

0.2a 

0.21a 

Argyropoulos  et  al.  [25] 

a  Estimated  value. 
b  Experimentally  measured  value. 


reported  in  the  literature  by  various  researchers.  Unfortunately, 
limited  experimental  data  are  available  for  thermal  conductivity 
and  contact  resistance  of  PEFC  material  to  support  the  math¬ 
ematical  models.  This  work  is  meant  to  provide  reliable  bench 
mark  data  for  modelers  to  accurately  predict  the  temperature  dis¬ 
tribution  in  a  fuel  cell.  In  this  paper,  an  experimental  technique 
based  on  a  steady  state  measurement  method  [19,26]  is  used 
to  accurately  measure  the  thermal  conductivity  of  dry  Nafion® 
membrane,  diffusion  media  (DM)  with  different  PTFE  content, 
and  from  various  manufacturers,  a  catalyst  layer,  and  the  contact 
resistance  between  the  diffusion  media  and  metal  plate  as  func¬ 
tion  of  temperature  and/or  compression  pressure.  Care  was  taken 
to  provide  accurate  assessment  of  the  experimental  uncertainty. 


Fig.  1.  Schematic  diagram  of  the  test  setup  to  measure  thermal  conductivity  and  thermal  contact  resistance  of  thin  films. 
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2.  Experimental  setup 

The  schematic  diagram  of  the  test  rig  for  the  thermal  contact 
resistance  and  through-plane  thermal  conductivity  measurement 
of  thin  films  is  shown  in  Fig.  1 .  In  this  configuration,  an  insulated 
long  circular  cross  section  of  standard  material  with  known  ther¬ 
mal  conductivity  is  uniformly  compressed  between  two  backing 
plates.  Both  backing  plates  are  maintained  at  a  different  constant 
temperature,  thus  acting  as  a  heat  source  and  sink,  to  generate 
the  heat  flow.  Thermocouple  arrays  are  inserted  in  the  cylindrical 
standard  material  above  and  below  the  test  sample.  The  test  sam¬ 
ple  with  unknown  thermal  conductivity  is  placed  between  the 
two  standard  material  cylinder  rods.  The  compression  pressure 
is  precisely  controlled  and  measured  with  a  load  cell. 

2.1.  Design  and  instrumentation 

The  test  setup  was  designed  based  on  two  criteria.  First  is  to 
achieve  nearly  one-dimensional  heat  transfer  in  the  axial  direc¬ 
tion.  And  second  was  to  achieve  a  temperature  drop  in  adjacent 
thermocouple  locations  that  is  comparable  in  magnitude  to  the 
temperature  drop  across  the  test  specimen,  so  that  the  signal  to 
noise  ratio  is  maximized.  Comprehensive  design  and  analysis 
were  performed  to  select  the  standard  material  based  on  its  ther¬ 
mal  conductivity,  to  determine  the  assembly  dimensions  and 
thickness  of  the  insulation.  Based  on  this  analysis,  aluminum 
bronze  (C613,  k  =  54.5 Wm-1  K-1)  was  chosen  as  the  stan¬ 
dard  material  above  and  below  the  specimen.  The  test  specimen 
was  placed  between  two  cylindrical  pieces  of  standard  mate¬ 
rial  7.62  cm  (3  in.)  in  length  and  5.08  cm  (2  in.)  in  diameter. 
Three  high  precision  thermocouples  (TMQSS-125-G-6,  Omega 
Engineering)  per  array  were  fixed  equidistantly  in  the  upper  and 
lower  pieces  of  the  standard  material.  The  temperatures  of  upper 
and  lower  compression  plates  were  monitored  by  a  T-type  fast 
response  self-adhesive  thermocouple  (SA1-T,  Omega  Engineer¬ 
ing).  A  load  cell  (LC304-1K,  Omega  Engineering)  was  placed 
between  the  two  bottom  plates  to  measure  the  compression 
pressure.  To  limit  radial  heat  leakage,  thermal  insulation  was 
provided  by  several  layers  of  polyethylene  foam  wrap  around  the 
standard  material/specimen.  The  steady  upper  and  lower  plate 
temperatures  were  maintained  by  passing  high  flow  rate  coolant 
fluid  through  channels  inside  the  plates  from  an  insulated  pro¬ 
grammable  temperature  controller  baths  (VWR  1 157P).  During 
experimentation,  the  entire  system  was  monitored  at  a  specified 
condition  for  8-10  h  to  ensure  a  true  steady  state.  Steady  state 
was  assumed  to  be  reached  when  the  temperature  fluctuation 
in  each  measured  point  in  the  system  was  less  than  0.5  °C  in 
30  min. 

2.2.  Measurement  and  uncertainty  analysis 

At  steady  state,  the  temperature  at  each  location  in  the  upper 
and  lower  thermocouple  arrays  was  measured  and  averaged  over 
a  time  span  of  20-30  min  for  the  analysis.  The  results  were  lin¬ 
early  extrapolated  to  the  test  specimen  edge,  to  determine  the 
temperature  drop  across  the  test  specimen.  Fig.  2  shows  a  plot  for 
an  experiment  with  127  p,m  (0.005  in.)  polyethylene  film  as  the 


Fig.  2.  Steady  state  temperature  profile  at  thermal  probe  location  for  validation 
experiment  with  127  p,m  (0.005  in.)  thick  polyethylene  sample  (vertical  line 
shows  the  location  of  thin  film). 


test  specimen.  With  the  knowledge  of  standard  material  thermal 
conductivity,  the  precise  heat  flux  across  the  specimen  can  be 
calculated  using  Fourier’s  law  [27]  in  the  axial  direction.  Since 
the  heat  flux  is  proportional  to  the  temperature  difference,  the 
through-plane  thermal  resistance  across  the  thin  film  test  spec¬ 
imen  can  be  easily  determined.  This  measurement  includes  the 
thermal  resistance  provided  by  test  specimen  and  thermal  con¬ 
tact  resistance  between  test  specimen  and  standard  material.  By 
performing  the  experiment  with  two  different  thickness  of  test 
specimen,  the  intrinsic  thermal  conductivity  (thermal  resistance 
of  specimen)  and  thermal  contact  resistance  between  the  speci¬ 
mens  can  be  determined  at  the  measured  compression  pressure 
and  temperature. 

The  thermal  contact  resistance  between  two  surfaces  is  a 
function  of  pressure  and  surface  roughness  [27].  The  test  sys¬ 
tem  allows  for  the  investigation  of  effect  of  pressure  on  thermal 
contact  resistance.  Surface  roughness  effects  were  minimized 
in  the  design  by  using  polished  metal  contact  surfaces.  A  dig¬ 
ital  display  connected  with  the  load  cell  allows  the  applied 
load  to  be  determined.  The  effect  of  temperature  on  the  intrin¬ 
sic  through-plane  thermal  conductivity  was  studied  by  varying 
the  heat  flow  through  the  test  specimen.  This  was  achieved  by 
altering  the  temperature  set  points  of  the  end  plates.  Temper¬ 
ature  controllers  can  provide  set  point  temperature  from  —15 
to  120  °C.  This  constrained  along-with  the  second  design  crite¬ 
ria  provides  a  minimum  and  maximum  average  temperature  of 
10-15  and  60-70  °C  for  the  test  specimen,  respectively. 

The  test  system  was  designed  to  minimize  the  uncertainty  in 
the  experiment.  Heat  loss  in  the  radial  direction  and  the  ther¬ 
mocouple  measurement  are  considered  to  be  the  major  potential 
sources  of  uncertainty.  Uncertainty  in  temperature  measurement 
is  accounted  for  probe  placement  in  the  hole  [28]  and  its  tol¬ 
erance.  Uncertainty  in  the  specimen  dimension  is  accounted 
for  the  manufacturing  tolerance.  A  detailed  uncertainty  analysis 
was  conducted  according  to  the  method  described  by  Kline  and 
McClintock  [29].  This  method  allows  estimation  of  the  uncer¬ 
tainty  of  experimental  measurements.  If  a  given  variable  G  is  a 
function  of  several  variables  so  that  G  =  G(x\ ,  xj,  ■  ■ .,  xn),  then 


M.  Khandelwal,  M.M.  Mench  /  Journal  of  Power  Sources  161  (2006)  1106-1115 


1109 


the  uncertainty  in  G,  providing  uncertainties  in  x\-xn  are  not 
coupled,  is  shown  as 


UG 


(1) 


where  uq  is  the  uncertainty  of  variable  G,  and  u\ -u„  are  the 
uncertainties  of  variable  x\ -xn ,  respectively.  For  the  thermal  con¬ 
ductivity  of  specimen,  all  dimensions  and  thermal  conductivity 
of  standard  material,  specimen  thickness,  and  temperature  mea¬ 
surement  were  considered  as  the  dependent  variables.  Based  on 
this  analysis,  the  estimated  uncertainty  in  the  thermal  resistance 
(thermal  conductivity)  of  the  specimens  should  be  less  than  15% 
of  measured  value,  with  most  uncertainty  as  a  result  of  radial  heat 
leak.  With  the  current  design  and  dimensions,  the  estimated  ratio 
of  radial  and  axial  flux  should  be  less  than  10%. 


3.  Result  and  discussion 

3.1.  Measurement  technique  validation 

The  test  system  and  measurement  technique  was  validated 
with  both  metal  and  thin  films.  Aluminum  bronze  (same  as 
the  standard  material)  samples  of  thickness  2.54  cm  (1  in.)  and 
5.08  cm  (2  in.)  were  used  to  validate  the  measurement  technique. 
The  thermal  conductivity  measured  was  8%  higher  than  the 
actual  value,  which  is  within  the  expected  experimental  uncer¬ 
tainty  and  is  likely  the  result  of  some  radial  heat  leak.  Polyethy¬ 
lene  films  of  thickness  127  p,m  (0.005  in.),  254  p,m  (0.010  in.) 
and  508  pun  (0.020  in.)  were  used  to  validate  the  test  setup  for 
thin  films.  The  calculated  thermal  conductivity  for  the  polyethy¬ 
lene  film  was  estimated  0.35  ±  0.05  W  m~ 1  K± 1 .  The  measured 
value  is  10%  higher  than  the  reported  value  in  the  literature 
0.32Wm_1K_l  [11,30],  and  lies  within  the  calculated  mea¬ 
surement  uncertainty.  An  overestimation  of  the  actual  thermal 
conductivity  is  expected,  due  to  the  radial  heat  loss.  However, 
since  the  exact  radial  heat  loss  varies  between  cases,  it  can¬ 
not  be  arbitrarily  adjusted  to  correct  in  the  thermal  conductivity 
value.  Nevertheless,  it  is  expected  that  all  values  reported  here 
are  5-10%  higher  than  the  actual  values  because  of  this  effect. 
Results  of  a  steady  state  temperature  profile  for  the  validation 


Fig.  3.  Measured  thermal  conductivity  as  a  function  of  temperature  for  dry 
Nation®  membrane  (using  Nation®  1 12  and  Nation®  1 17),  error  bar  calculated 
at  15%. 


experiment  for  a  127  pm  (0.005  in.)  thin  polyethylene  sample 
are  shown  in  Fig.  2.  The  experimentally  measured  thermal  con¬ 
ductivities  of  various  commonly  used  fuel  cell  components  is 
summarized  in  Table  2,  and  discussed  in  the  following  sections. 

Repeatability  of  measurement  is  estimated  to  be  15-20% 
based  on  repeated  test  on  Nafion®  membrane  at  high  compres¬ 
sion  pressure  (~2MPa).  However,  additional  variation  can  be 
expected  at  reduced  compression  pressure  due  to  non-uniform 
contact. 

3.2.  Nafion®  membrane 

The  thermal  conductivity  of  DuPont  Nafion®  membrane 
(provided  by  IonPower  Inc.)  over  a  range  of  temperatures  from 
17  to  65  °C  is  plotted  in  Fig.  3.  To  estimate  dry  Nafion®  ther¬ 
mal  conductivity  (at  atmospheric  condition),  Nafion®  112  and 
Nafion®  117  samples  were  tested  as  delivered,  with  no  exter¬ 
nal  humidification.  At  room  temperature,  dry  Nafion®  thermal 
conductivity  was  found  to  be  0.16  ±  0.03  W  nr-1  K-1.  At  the 
highest  measured  temperature  of  65  °C,  the  thermal  conductiv¬ 
ity  of  Nafion®  was  determined  to  be  0.13  ±0.02  W m_1  K-1. 
At  80  °C,  the  thermal  conductivity  can  be  linearly  extrapolated 
as  0. 1 1  ±  0.02  W  m_  1  K~ 1 .  The  measured  thermal  conductivity 


Table  2 

Experimentally  measured  thermal  conductivity  of  PEFC  components 


Material 

Measured  thermal 
conductivity,  k  (W  m-1  K-1) 

Reported/estimated  value 
in  literature  (W  m-1  K-1) 

DuPont  Nafion®  membrane  (at  30  °C) 

0.16  ±  0.03 

0.1-0.95 

Toray  carbon  fiber  paper  diffusion  media  (TGP-H-60  at  26  °C) 

1.80  ±  0.27 

1.7  [37] 

SIGRACET®  0  wt.%  PTFE  carbon-fiber  paper  diffusion  media  (AA  series  at  56  °C) 

0.48  ±  0.09 

Not  available 

SIGRACET®  5  wt.%  PTFE  carbon-fiber  paper  diffusion  media  (BA  series  at  58  °C) 

0.31  ±  0.06 

Not  available 

SIGRACET®  20  wt.%  PTFE  carbon-fiber  paper  diffusion  media  (DA  series  58  °C) 

0.22  ±  0.04 

Not  available 

E-Tek  ELAT®  diffusion  media  (LT1200-W  at  33  °C) 

0.22a  ±  0.04 

Not  available 

Catalyst  layer  (0.5  mg  cm-2  platinum  on  carbon) 

0.27b  ±  0.05 

Not  available 

a  Thermal  conductivity  based  on  single  and  double  layer  of  diffusion  media  as  test  specimen  (neglecting  contact  resistance  between  the  two  layers  of  diffusion 
media). 

b  Effective  thermal  conductivity  (includes  thermal  contact  resistance  with  diffusion  media). 


1110 


M.  Khandelwal,  M.M.  Mench  /  Journal  of  Power  Sources  161  (2006)  1106-1115 


is  in  good  agreement  with  the  experimental  observation  of  Vie 
and  Kjelstrup  [22],  and  Burford  and  Mench  [23]  as  shown  in 
Table  1. 

It  was  observed  that  the  thermal  conductivity  of  Nafion® 
decreases  with  increasing  temperature,  where  the  test  tempera¬ 
ture  is  defined  as  the  average  temperature  across  the  test  spec¬ 
imen.  Although  the  downward  trend  in  thermal  conductivity  is 
nearly  within  the  error  bar  limit,  this  behavior  of  Nafion®  may  be 
explained  on  the  basis  of  its  structure  and  morphology.  Kenneth 
and  Moore  [31],  and  other  researchers  [32,33]  have  presented  a 
comprehensive  review  on  the  fundamental  structure  and  prop¬ 
erties  of  Nafion®  perfluoronate  materials  and  its  applications. 
Using  small  angle  X-ray  scattering  (SAXS)  and  wide  angle  X- 
ray  diffraction  (WAXD),  Gierke  and  co-workers  [34,35]  also 
examined  and  compared  morphological  features  of  Nafion®  and 
attributed  these  features  to  a  crystalline  structure  within  the  flu¬ 
orocarbon  matrix.  For  a  crystalline  polymer,  Choy  et  al.  [  1 1,36] 
has  explained  an  inverse  relation  between  temperature  and  ther¬ 
mal  conductivity  based  on  the  phonon  transport.  They  reported 
that  the  thermal  conductivity,  k ,  of  an  insulating  material  is  given 
as  [11,18]: 

k  =  '  cvl  (2) 

3 

where  c  is  the  phonon  heat  capacity,  v  the  phonon  velocity  and  /  is 
the  phonon  mean  free  path.  The  collision  probability  of  phonons 
increases  at  high  temperature  due  to  an  increase  in  the  number 
of  participating  phonons.  Consequently,  the  mean  free  path  / 
decreases,  while  the  c  and  v  remain  almost  constant.  Thus,  the 
thermal  conductivity,  which  is  proportional  to  the  mean  free 
path,  decreases  with  an  increase  in  the  temperature  for  a  crys¬ 
talline  insulating  material,  such  as  Nafion®. 

The  effect  of  humidity  on  the  thermal  conductivity  of  the 
Nafion®  membrane  can  be  theoretically  estimated  by  assuming 
the  humidified  membrane  as  mixture  of  water  and  dry  mem¬ 
brane.  Using  volume  averaging,  the  thermal  conductivity  of  a 
humidified  membrane  can  be  estimated  as 

T  (hy//Pv/)kw  ±  (l/csch~  Mw)knf 

k  =  - - -  (3) 

±  Aw/pw 


where  Xw  is  the  water  content  in  Nafion®  (mol  HbO/mol  SO3- 1 ), 
pw  the  density  of  water  (kgm-3),  kw  the  thermal  conduc¬ 
tivity  of  water  (Wm-IK-1),  Mw  the  molecular  weight  of 
water  (kgkmol-1),  cS03-  the  molar  concentration  of  SO3-  ion 
(moll-1)  and  knf  is  the  thermal  conductivity  of  dry  Nafion® 
membrane  (W  m-1  K-1).  Note  that  the  non-homogenous  water 
distribution  in  the  membrane  may  cause  error  in  this  estimation. 
Nevertheless,  it  serves  to  bound  the  expected  thermal  conduc¬ 
tivity  of  a  humidified  membrane  for  calculation  and  comparison 
purposes.  From  Eq.  (3),  the  estimated  variation  of  thermal  con¬ 
ductivity  of  humidified  Nafion®  membrane  with  temperature 
was  calculated,  and  is  shown  in  Fig.  4  for  different  humidity 
ratios.  For  a  fully  humidified  membrane,  the  thermal  conductiv¬ 
ity  was  estimated  to  be  0.29  ±  0.03  W  m-1  K-1  at  65  °C.  It  is 
important  to  observe  that  water  thermal  conductivity  increases 
with  temperature;  however  the  inverse  relation  is  true  for  dry 
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Fig.  4.  Estimated  thermal  conductivity  of  Nafion®  1 100  EW  at  different  humid- 
ity  ratios.  The  thermal  conductivity  variation  of  pure  water  is  shown  as  an  upper 
bound  for  the  theoretical  moist  Nafion®  thermal  conductivity. 


Nafion®.  Hence,  for  the  case  of  80%  or  fully  humidified  con¬ 
dition  (100%),  the  membrane  thermal  conductivity  is  uniform 
with  temperature  within  the  experiment  tolerance.  The  thermal 
conductivity  of  pure  liquid  water  is  shown  in  Fig.  4  as  an  upper 
bound  on  the  maximum  possible  thermal  conductivity. 

3.3.  Diffusion  media 

Toray  carbon  paper  “TGP-H”  [37]  (0wt.%  PTFE)  diffu¬ 
sion  media  were  investigated  to  measure  the  thermal  con¬ 
ductivity  and  thermal  contact  resistance  with  the  aluminum 
bronze  standard  material  under  various  temperature  and  com¬ 
pression  pressures.  Two  diffusion  media  samples  with  thickness 
of  190  |xm  (TGP-H-60,  78%  porosity)  and  280  |xm  (TGP-H- 
90,  78%  porosity)  were  used  in  the  experiment.  The  measured 
intrinsic  material  thermal  conductivity  at  26  °C  was  found  to 
be  1.80  ±  0.27  Wm-1  K-1,  which  is  in  close  agreement  with 
the  value  provided  in  the  specification  of  the  material  from  the 
manufacturer  of  1 .70  W  m- 1  K- 1  [37] .  Variation  of  thermal  con¬ 
ductivity  of  Toray  carbon  paper  (TGP-H-60)  with  temperature 
is  shown  in  Fig.  5  and  it  decreases  from  1.80  ±  0.27  W  m-1  K-1 
at  26  °C  to  1.24±  0.19  W  m-1  K-1  at73°C.  This  inverse  rela¬ 
tion  of  thermal  conductivity  with  temperature  for  carbon  paper 
may  be  explained  based  on  the  presence  of  carbonized  thermo 
setting  resin  [41],  which  acts  as  a  binder  in  the  diffusion  media. 
Based  on  the  phonon  transport  [11,18],  thermal  conductivity  of 
these  thermo-setting  polymer  decreases  with  increasing  temper¬ 
ature.  For  diffusion  media  without  any  binder  material,  insignif¬ 
icant  changes  (±0.04  W  m- 1  K- 1 )  in  thermal  conductivity  were 
observed  for  carbon  cloth  (ELAT®  LT  1200-W)  in  the  temper¬ 
ature  range  of  30-60  °C. 

Variation  of  thermal  contact  resistance  of  Toray  carbon  paper 
with  aluminum  bronze  material  for  different  compression  pres¬ 
sure  was  calculated  and  is  shown  in  Fig.  6.  In  this  experiment, 
the  initial  pressure  was  gradually  increased  from  0.4  to  2.2  MPa, 
and  then  reduced  to  0.4  MPa.  Hysteresis  is  observed  due  to 
diffusion  media  fiber  compression  and  deformation.  An  asymp- 
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Fig.  5.  Measured  thermal  conductivity  as  a  function  of  temperature  for  Toray 
carbon  paper  (TGP-H-60),  error  bar  calculated  at  15%. 


totic  approach  to  a  minimum  value,  corresponding  to  a  fully 
compressed  diffusion  media  surface,  is  seen  near  2  MPa  pres¬ 
sure  ranges,  which  corresponds  to  a  thermal  contact  resistance 
of  2.0  x  10~4  m2  KW"1.  Contact  resistance  is  also  function  of 
type  of  the  material.  So,  the  reported  value  here  may  differ  from 
the  contact  resistance  with  graphite  or  other  cell  materials  with 
different  surface  finish.  The  point  is  to  show  the  hysteresis  that 
exists  as  a  function  of  compression  pressure. 

3.4.  Effect  ofPTFE  content  on  diffusion  media  thermal 
conductivity 

To  study  the  effect  of  hydrophobic  additive  polytetrafluo- 
roethylene  (PTFE)  on  the  diffusion  media  thermal  conductiv¬ 
ity,  experiments  were  performed  on  SIGRACET®  AA  (Owt.% 
PTFE),  BA  (5  wt.%  PTFE)  and  DA  (20  wt.%  PTFE)  series  dif¬ 
fusion  media  [38],  with  thicknesses  of  190  and  280  p.m.  The 
variation  of  thermal  conductivity  with  PTFE  content  for  these 
SIGRACET®  diffusion  media  is  shown  in  Fig.  7.  It  was  observed 


Fig.  6.  Variation  of  thermal  contact  resistance  between  Toray  carbon  paper 
and  smooth  aluminum  bronze  with  compression  pressure,  TGP-H-60  (0.19  mm 
thick)  and  TGP-H-90  (0.28  mm  thick)  with  78%  porosity. 


with  wt.%  PTFE  content  at  2  MPa  compression  pressure  and  56-58  °C. 

that  thermal  conductivity  decreases  with  increasing  PTFE  con¬ 
tent.  For  low  PTFE  content  (5  wt.%),  the  thermal  conductivity 
(0.31  ±  0.06  Wm-1  K-1)  is  reduced  by  35%  as  compared  to  the 
untreated  diffusion  media  (0.48  ±  0.09  Wm-1  K-1).  However, 
with  increase  in  weight  percent  of  PTFE  content  from  5  to  20, 
the  diffusion  media  thermal  conductivity  was  further  reduced 
only  by  29%  (0.22  ±  0.04  W  m“ 1  K“ 1 ).  Therefore,  with  further 
increase  in  the  weight  percent  of  PTFE  the  thermal  conductivity 
of  PTFE  treated  diffusion  media  could  approach  an  asymptotic 
value.  The  porosity  value  of  all  SIGRACET®  diffusion  media 
was  82-85%,  implying  some  of  the  high  thermal  conductive 
carbon  fibers  (£carbon  =  129  W  m-1  K-1  [39])  are  coated  by  low 
thermally  conducting  PTFE  material  (£ptfe  =  1 1 .7  W  m_  1  K“ 1 
[40]).  This  could  be  attributed  as  one  of  the  possible  reasons 
for  reduction  in  thermal  conductivity  of  PTFE  treated  diffusion 
media.  Scanning  electron  microscopy  (SEM)  images  of  PTFE 
treated  and  untreated  SIGRACET®  diffusion  media  are  given 
in  Fig.  8.  Other  work  [41]  has  also  identified  the  PTFE  distri¬ 
bution  preferentially  near  the  diffusion  media  surface.  Due  to 
the  coverage  of  these  carbon  fibers,  the  thermal  resistance  is 
increased,  which  leads  to  the  reduction  in  the  effective  thermal 
conductivity  of  PTFE  treated  diffusion  media.  In  theory,  PTFE 
addition  can  actually  increase  the  thermal  conductivity  of  the 
DM  if  the  PTFE  displaces  a  significant  amount  of  lower  ther¬ 
mal  conductive  air.  However,  this  was  not  the  case.  The  thermal 
contact  resistance  was  determined  to  be  1.3  x  10  ,  1.5  x  10~4 

and  1.2  x  10“4m2KW“1  for  SIGRACET®  untreated,  5  and 
20  wt.%  PTFE  content  diffusion  media,  respectively  at  2.2  MPa. 

3.5.  Catalyst  layer 

The  catalyst  layer  which  is  a  mixture  of  electrolyte  and 
0.5  mg  cm-2  platinum  powder  loading  supported  by  large  car¬ 
bon  particles,  was  hot  pressed  on  both  sides  of  a  polymer  elec¬ 
trolyte  membrane.  This  three-layer  membrane  structure  is  often 
referred  as  membrane  electrode  assembly  (MEA).  To  determine 
the  through-plane  thermal  resistance  of  the  catalyst  layer,  exper¬ 
iments  were  performed  on  a  MEA  (Nafion®  112  membrane. 
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Fig.  8.  Scanning  electron  microscopy  (SEM)  image  of  SIGRACET®  PTFE  treated  and  untreated  diffusion  media  at  different  magnification:  (a)-(b)  untreated  diffusion 
media  (type  24  AA);  (c)-(d)  20  wt.%  PTFE  diffusion  media  (type  24  DA). 


cathode/anode  catalyst  layer:  25  |xm,  manufacturer:  Lynntech, 
provided  by  FuelCellStore.com)  having  a  diffusion  media  (Toray 
carbon  paper:  TGP-H-60)  on  both  side.  From  previous  exper¬ 
iments,  the  membrane  thermal  conductivity,  diffusion  media 
thermal  conductivity  and  its  contact  resistance  with  the  standard 
material  are  known.  So,  the  combined  thermal  resistance  of  the 
catalyst  layer  and  contact  resistance  with  the  diffusion  media 
can  be  determined.  This  value  was  termed  the  effective  thermal 
resistance  of  the  catalyst  layer.  Variation  of  the  effective  thermal 
contact  resistance  with  compression  pressure  for  a  catalyst  layer 
is  shown  in  Fig.  9.  An  asymptotic  region  is  approached  at  com- 


Fig.  9.  Variation  of  effective  thermal  resistance  (including  thermal  contact  resis¬ 
tance  with  Toray  carbon  paper)  of  0.5  mg  cm-’  platinum  loaded  catalyst  layer 
with  compression  pressure. 


pression  pressure  ~2  MPa  which  corresponds  to  thermal  contact 
resistance  1  x  10-4m2KW-1.  This  thermal  resistance  value 
may  be  very  helpful  for  non-isothermal  mathematical  models 
which  presently  ignore  the  catalyst  layer  contribution  to  the  ther¬ 
mal  resistance.  The  equivalent  effective  thermal  conductivity  of 
the  catalyst  layer  was  also  estimated  as  0.27  ±  0.05  W  m-1  K-1, 
which  is  similar  in  magnitude  to  that  of  the  diffusion  media. 

3.6.  Theoretical  study  of  temperature  drop  in  a  fuel  cell 

A  simple  analytical  model  based  on  one  dimensional  con¬ 
duction  heat  transfer  was  developed  to  predict  the  typical  tem¬ 
perature  drop  in  the  fuel  cell  components  using  the  new  thermal 
transport  parameter.  This  one-dimensional  approach  is  reason¬ 
able  as  a  first  approximation  to  investigate  the  effect  of  material 
properties  on  the  temperature  distribution  in  a  fuel  cell. 

The  simplified  ID  model  accounts  for  the  heat  transfer  from 
the  membrane,  anode/cathode  diffusion  media,  catalyst  layer 
and  backing  plate,  as  shown  in  Fig.  10.  Steady  state  is  assumed 
and  conduction  heat  transfer  from  the  flow  channel  is  neglected, 
since  it  is  typically  relatively  small.  The  heat  source  is  consid¬ 
ered  from  Joule  heating  in  the  membrane.  In  the  catalyst  layers, 
the  heat  source  is  calculated  from  the  respective  half-cell  reac¬ 
tion  entropic  loss,  activation,  concentration  over-potential,  and 
joule  heating.  The  governing  equation  in  each  component  can 
be  written  from  the  Fourier’s  law  [27]: 


where  k  is  the  thermal  conductivity  of  the  material  ( W  m  1  K  1 ), 
T  the  temperature  (K)  and  q'"  (W  m-3)  is  the  volumetric  heat 
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Fig.  10.  Schematic  of  polymer  electrolyte  fuel  cell  for  the  simplified  ID  ana¬ 
lytical  conduction  model  for  thermal  transport. 


Fig.  1 1 .  Estimated  temperature  drop  in  fuel  cell  component  for  current  density 
i  =  0. 1  and  1 .0  A  cm-2  for  Toray  carbon  paper  and  SIGRACET®  5%  PTFE  as  the 

diffusion  media.  ( — )  i  =  0.1  A  cm-2  (Toray),  ( - )  i=  1.0  A  cm-2  (Toray), 

( - )  i  =  0. 1 A  cm-2  (SIGRACET®);  ( - )i  =  1.0  A  cm"2  (SIGRACET®). 


source.  A  constant  temperature  boundary  condition  is  specified 
on  both  the  backing  plate,  and  temperature  and  heat  flux  conti¬ 
nuity  is  imposed  on  each  interface  in  the  domain.  Thermal  con¬ 
tact  resistance  is  also  accounted  on  the  backing  plate-diffusion 
media  interface  and  the  value  is  used  from  our  experimental 
result.  The  heat  generation  in  catalyst  layer  can  be  written  as 
[42]: 


q 


m 


i 

^cl 


(5) 


where  n  is  the  number  of  electrons  per  mole  of  the  reactants, 
AS  the  half-cell  change  in  the  entropy  (Jmol-1K-1),  i  the 
current  density  (A  cm-2),  ?)  the  voltage  drop  (over-potential) 
due  to  activation,  concentration  and  ohmic  resistance  (V),  fcj 
the  thickness  of  the  catalyst  layer  (|xm)  and  F  is  the  Faraday 
constant.  The  entropy  change  at  standard  state  with  platinum 
catalyst  is  taken  as  AS  =  0.104Jmol-1  K  for  the  anode  side, 
and  A S=  —326.36  J  mol-1  K-1  for  the  cathode  side  [43].  The 
activation  over-potential  was  calculated  based  on  typical  Tafel 
kinetics  for  a  Pt-electrode  at  the  normal  operating  condition.  For 
the  membrane  ohmic  resistance  (Joule  heating)  is  accounted  for 
volumetric  heat  generation. 

Based  on  Eqs.  (4)  and  (5),  a  temperature  profile  in  the  fuel 
cell  components  was  obtained  and  shown  in  Fig.  1 1  for  two  dif¬ 
ferent  current  densities  of  0.1  and  1  A  cm-2  with  SIGRACET® 
and  Toray  carbon  paper  diffusion  media.  The  model  parameters 
are  listed  in  Table  3.  For  a  200  |xm  thick  diffusion  media,  5  wt.% 
PTFE  content  SIGRACET®  has  an  almost  3^4  °  C  drop  in  the  dif¬ 
fusion  media  as  compared  to  1-2  °C  temperature  drop  in  Toray 
carbon  paper  for  i=  1  A  cm-2.  However,  temperature  is  almost 
uniform  for  lower  current  density,  i.e.  /  =  0. 1  A  cm-2.  From  the 
scale  analysis,  7-8  °C  temperature  drop  can  be  expected  for 
400  p,m  thick  diffusion  media  at  i  =  1.0  A  cm-2.  A  temperature 
drop  of  10-12  °C  was  also  measured  by  Burford  and  Mench 
[23],  and  estimated  by  Hwang  [3]  for  much  thicker  diffusion 


Table  3 

Geometric  and  material  parameters  used  in  simplified  analytical  model 


Parameter 

Value 

Membrane  thickness 

50  jjim 

Anode/cathode  diffusion  media  thickness 

200  |jim 

Anode/cathode  catalyst  layer  thickness 

25  p,m 

Backing  plate  thickness 

2  mm 

Membrane  thermal  conductivity 

0. 12  W  m_I  K"1 

Membrane  ionic  conductivity 

8  3  m-1  [2,4] 

Anode/cathode  diffusion  media  thermal 

1.7  Wm^1  KT1  (Toray), 

conductivity 

0.31  Wm-1  KT1  (SIGRACET® 

Anode/cathode  catalyst  layer  thermal 

5  wt.%) 

0.27  Wn-T1  KT1 

conductivity 

Backing  plate  (BP)  thermal  conductivity 

55Wm-‘  KT1 

Thermal  contact  resistance  BP-diffusion 

2  x  10-4  m2  KW_1 

media 

Backing  plate  end  temperature 

80  °C 

media,  and  is  consistent  with  the  current  model  results.  This 
temperature  drop  has  a  significant  affect  on  the  water  saturation 
curve. 

All  experiments  performed  in  the  current  work  were  to  mea¬ 
sure  thermal  conductivity  in  the  through-plane  direction  for  fuel 
cell  components.  Due  to  the  highly  anisotropic  nature  of  porous 
diffusion  media,  thermal  conductivity  in  the  in-plane  direction 
also  plays  a  very  crucial  role  in  the  heat  transfer  in  PEFC  and 
will  be  investigated  in  future  study. 

4.  Conclusion 

An  experimental  technique,  based  on  the  steady  state  method 
was  used  to  precisely  determine  the  though-plane  thermal  con¬ 
ductivity  and  contact  resistance  of  varying  fuel  cell  components, 
including  the  electrolyte,  catalyst  layer  and  several  types  of  dif- 
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fusion  media.  The  following  conclusions  can  be  made  regarding 
this  work: 

(i)  Thermal  conductivity  of  dry  Nafion®  membrane  was 
measured  to  be  0.16  ±  0.03  W  m_1  K-1  at  30  °C,  and 
the  extrapolated  thermal  conductivity  at  80  °C  was 
0.11  ±  0.02  W  m-1  K  ,  decreases  with  increase  in  tem¬ 
perature. 

(ii)  The  effect  of  humidity  on  thermal  conductivity  of  the 
Nafion®  membrane  was  estimated  using  a  volume  average 
technique.  The  thermal  conductivity  of  Nafion®  increases 
with  the  water  content.  For  a  fully  humidified  mem¬ 
brane,  the  thermal  conductivity  was  estimated  to  be 
0.29  ±  0.03  W  m~ 1  K“ 1  at65°C. 

(iii)  The  experimentally  measured  thermal  conductivity  for 
Toray  carbon  paper  was  1.80  ±  0.27  W  m-1  K-1  at  26  °C 
and  it  decreases  with  increase  in  temperature.  The  tightly 
compressed  thermal  contact  resistance  for  Toray  carbon 
paper  with  smooth  aluminium  bronze  material  was  esti¬ 
mated  as  2.0  x  10-4  m2  KW"1  at  a  compression  pressure 
of  2  MPa  and  over. 

(iv)  The  thermal  conductivity  of  PTFE  treated  diffusion  media 
decreases  with  increase  in  the  PTFE  content.  The  thermal 
conductivity  of  S1GRACET®  untreated  (0wt.%),  5  and 
20wt.%  PTFE  content  diffusion  media  was  measured  to 
be  0.48  ±0.09,  0.31  ±0.06  and  0.22  ±  0.04  Wm-1  KT1, 
respectively. 

(v)  The  thermal  contact  resistance  of  SIGRACET®  diffu¬ 
sion  media  with  aluminium  bronze  material  is  almost 
constant  with  varying  PTFE  content  and  is  order  of 
1.2-1. 5  x  10-4  m2  KW-1  at  2.2  MPa. 

(vi)  The  effective  thermal  conductivity  of  a  catalyst  layer  hav¬ 
ing  0.5  mg  cm-2  platinum  on  carbon  was  estimated  to  be 
0.27  ±  0.05  Wm-’K"1 

(vii)  A  simple  one-dimensional  analytical  model  was  used  to 
predict  the  temperature  drop  of  3-4  °C  was  estimated 
for  200  p,m  thick  SIGRACET®  5wt.%  diffusion  media 
at  1  Acm~2  80  °C  operating  condition. 

This  study  should  be  useful  to  resolve  discrepancies  in  miss¬ 
ing  values  of  thermal  conductivity  in  the  literature  of  various 
fuel  cell  materials  and  help  modelers  to  estimate  the  tempera¬ 
ture  distribution  in  an  operating  PEFC. 
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